Coupled opto-electronic simulation of organic bulk-heterojunction solar cells: 
parameter extraction and sensitivity analysis 
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A comprehensive opto-electronic device model for organic bulk-heterojunction solar cells is pre- 
sented. First the optical in-coupling into a multilayer stack is calculated. From the photon ab- 
sorption profile a charge transfer (CT) exciton profile is derived. In this study we consider the 
Onsager-Braun mechanism to calculate the dissociation of the CT-excitons into free charge carri- 
ers. These free charge carriers then migrate towards the electrodes under the influence of drift and 
diffusion. A general problem arising in computer simulations is the number of material and device 
parameters, which have to be determined by dedicated experiments and simulation-based parameter 
extraction. In this study we analyze measurements of the short-circuit current dependence on the 
active layer thickness and current- voltage curves in poly(3-hexylthiophene):[6,6]-phenyl-Cei-butyric 
acid methyl ester (P3HT:PCBM) based solar cells. We have identified a set of parameter values 
including dissociation parameters that describe the experimental data. The overall agreement of our 
model with experiment is good, however a discrepancy in the thickness dependence of the current- 
voltage curve questions the influence of the electric field in the dissociation process. In addition 
transient simulations are analyzed which show that a measurement of the turn-off photocurrent can 
be useful for estimating charge carrier mobilities. 

Keywords: organic photovoltaic, organic solar cell, P3HT:PCBM, numerical simulation, sensitivity, param- 
eter extraction, mobility measurement 
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I. INTRODUCTION 



After the huge progress organic light-emitting diodes 
have made in the last years, organic photovoltaic (OPV) 
devices attract more and more interest. State of the 
art OPV devices yield an energy conversion efficiency of 
around 6% to 7% for single junction cells El as well as 
tandem cells p|. This is much less compared to other 
already established photovoltaic power conversion tech- 
niques which have an efficiency above 10% or even above 
20% for crystalline silicon cells. Nevertheless OPV de- 
vices have several advantages like the possibility of low- 
cost production room temperature processing and 
thin film structures. The latter two make it possible to 
use flexible substrates and thus the production of flexi- 
ble solar cells. This leads to low-cost, flexible and trans- 
portable energy generators. Our research focuses on the 
influence of optical and electrical parameters on device 
performance. 

The field of OPV 's can be separated into the planar- 
heterojunction devices where the donor and acceptor ma- 
terials are deposited one after the other, mostly by evap- 
oration of small molecules and the bulk-heterojunction 
(BHJ) devices, where the two organic materials are di- 
luted in the same solvent and spin coated as one layer. 
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The advantage of the BHJ structure is that most of the 
generated excitons reach a nearby donor-acceptor inter- 
face where they dissociate into free charge carriers. This 
efficient exciton harvesting leads to higher power conver- 
sion efficiencies for BHJ devices. 

Numerical models for organic BHJ devices provide in- 
sight into operating mechanisms and allow for device 
structure optimization. The simulation of OPV devices 
can be separated into two parts, firstly there is the in- 
coupling of light into a multilayer structure and secondly 
the extraction of charges which needs an electrical model. 
The absorption of light within the multilayer structure is 
a crucial process and thus one of the main areas of numer- 
ical simulations. Harrison et al. [4[ used analytical mod- 
els and compared them to measured photocurrent action 
spectra. Later a transfer matrix formalism has been used 
by Petterson et al. [f| to calculate the absorbed optical 
energy within the multilayer structure. This formalism is 
widely used to optimize layer structures for single junc- 
tion [y] as well as tandem cells 0] ■ Though optical mod- 
els are very reliable they are only able to give an upper 
limit of the achievable photocurrent. Therefore electri- 
cal models have been developed to study charge genera- 
tion and transport losses in planar heterojunction [8| as 
well as BHJ devices [I EI El El . These models differ 
in their choice of ingredients, definition of the boundary 
conditions as well as in their methods to solve the drift- 
diffusion equations. Coupled opto-electronic models have 
been introduced by Kirchartz et al. EIIj Lacic et al. EH 
and Kotlarski et al. EH' where they connected a drift- 
diffusion solver with an optical thin-film simulator. All 
the electrical models mentioned are ID-models and they 
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FIG. 1: Electromagnetic field penetration plot (top) is calcu- 
lated using a transfer matrix formalism. This field penetra- 
tion is then used to derive the photon absorption rate profile 
(bottom). 

treat every layer as an effective medium, i.e. a homoge- 
nous layer. 

In this study we present a comprehensive numerical 
model which considers optical as well as electrical effects. 
This device model has originally been developed for the 
analysis of organic-light emitting diodes [l6l . [l7l . [hSj and 
was commercialised afterwards. This software, SETFOS 
[l9| , has been extended for the simulation of organic BHJ 
solar cells. 



II. DESCRIPTION OF THE NUMERICAL 
DEVICE MODEL 

A. Optical modeling 

Organic photovoltaic cells consist of a multilayer struc- 
ture, with layer thicknesses on the order of the wave- 
length of the incident light, which is shorter than the 
coherence length of sunlight. Therefore an OPV is an 
optically coherent device. This gives rise to interference 
effects inside the device which can be exploited to in- 
crease the maximum photocurrent by carefully tuning 
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FIG. 2: Reflectance and fractional absorbance for each layer 
in the solar cell. 

the thicknesses of the individual layers [2(| [E[ [ioj . An 
approach to find the optimal layer thicknesses from an 
optical point of view is to calculate the electromagnetic 
field inside the multilayer stack using a transfer matrix 
formalism. This approach has been used in the optical 
simulation of OLED's [3 QJ] and OPV's @ and is well 
established in the field of organic electronics. With this 
formalism the normalized optical electric field E norm is 
calculated throughout the multilayer stack (cf. Fig. [TJ. 
For plane waves, the optical field intensity within the 
stack is given by 



I=-e Q cn\E\ 2 . (1) 

In the transfer matrix formalism we normalize the in- 
cident electric field such that the incoming electric field 
amplitude is 1. Therefore the electric field calculation 
gives the field strength relative to the field strength in 
the surrounding medium E . This is used to calculate 
the absolute energy flux by looking at the ratio 

I_ = ±e cn\E\* 
Iq ±e cn \E \ 2 

of the illumination intensity Io in the surrounding 
medium and device-internally /. In this study the AM 
1.5 spectrum is used for Iq. The calculation of the device- 
internal intensity 

I — — \E norm \ 2 Io, (3) 
n 

is done wavelength and position dependent, from which 
the density of absorbed photons per second can be de- 
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FIG. 3: Relevant processes for CT-exciton modeling. 



rived via 



_ alX 

^photons — ~~ T i 



(4) 



where a stands for absorption coefficient which is given 
by a = and k stands for the complex part of the re- 
fractive index (i.e. the extinction coefficient). After an 
integration over the illumination spectrum one ends up 
with a photon absorption rate profile (cf. Fig. [1] (bot- 
tom)). The integration is done in the absorbing spectral 
range from 380 to 780 nm. The absorbance of the active 
layer, which is the ratio of absorbed to the illumination 
light intensity, drops to zero above 650 nm as is shown 
in Fig. [5J It is obvious that the thickness and refractive 
index is required for each layer for accurate simulation 
results. The n and k data were taken from Monestier et 
al. M. 



B. Electrical modeling 

To simulate an OPV device electrically, several pro- 
cesses have to be taken into account. First charge- 
transfcr-exciton generation and dissociation, then charge 
drift and diffusion, and at last the charge extraction at 
the electrodes need to be considered. These processes 
will be discussed in this section. 



1. Charge-transfer- exciton dissociation 

The absorption of a photon generates an exciton which 
then diffuses to a donor-acceptor interface. At this in- 
terface the exciton dissociates and generates a charge- 
transfer-exciton (CT-excitons). In this model we assume 
that a photon directly generates a CT-exciton. Therefore 
only a continuity equation for CT-excitons is needed (cf. 
Fig. 

dS{x) 



dt 

r e ffr(x)p(x)n(x) - kfS(x) - k d S(x) + g e ffG op t(x). 



CT-exciton density S is fed by two processes: Firstly 
an absorbed photon directly generates a CT-exciton. 
Gopt(x) stands for the photon absorption profile shown 
in figure [1] (bottom), whereas ff e // is a photon-to-CT- 
exciton conversion efficiency. We set g e // to zero out- 
side of the active layer, which is P3HT:PCBM in the 
device under consideration here. Secondly recombina- 
tion of free charge carrier pairs generates a CT-exciton. 
This recombination process is treated by the Langevin 
theory which will be discussed in the next section when 
the drift-diffusion model is explained. There are also two 
CT-exciton state depletion mechanisms kf and ka which 
stand for the decay of a CT-state and the dissociation 
respectively as described by the Onsager-Braun theory. 
This continuity equation for CT-excitons is coupled with 
the drift-diffusion model which will be discussed later. 
CT-exciton dissociation has been studied by L. Onsager 
pH | [22I ]. The initial pair separation distance is overes- 
timated when measurements are compared with the On- 
sager theory. Therefore C. L. Braun [23[ extended the 
theory by Onsager. In this model a finite decay rate of 
the CT-exciton has been considered. The dissociation 
probability P is then given by 



P(T,E,k f ,a) 



k d (T,E,a) 



A£{a) 



b(T,E) 



kd + kf 
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J x [2y/-2b(T,E) 
y/-2b(T,E) 



(5) 



(6) 

In equation ([5]) the mobility is the sum of the elec- 
tron and hole mobility fj, = ji e + fXh, ^i(x) is the Bessel 
function of the first kind of order 1, a is the initial pair 
separation distance of the CT-exciton, the pair binding 
energy A£(a) is calculated under the assumption that 
CT-excitons have the same dependence of the binding 
energy on the separation distance as ion pairs. E stands 
for the electrical field which is calculated position depen- 
dent using the drift-diffusion model explained in the next 
section. 

To find an analytical model, the Onsager-Braun theory 
makes some simplifications. It does not consider disorder 
in energy levels as well as in the spatial distribution of the 
hopping sites which is present in an amorphous organic 
semiconductor. Koster et al. Q and Pan et al. [24[ in- 
troduced spatial disorder by integrating the dissociation 
probability over distribution of pair separations a. A 
further simplification of the Onsager-Braun model is the 
assumption that the dissociation takes place in a homo- 
geneous material, neglecting the presence of donors and 
acceptors. Veldman et al. [25| proposed that the mobility 
used in the Onsager-Braun model is higher than the effec- 
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tive mobility measured in the bulk. They attributed this 
to local crystallization of the PCBM molecules. This ex- 
tension allows to model CT-exciton lifetimes in the order 
of 1CP 8 s still having a dissociation efficiency over 90%. 
This extension is not considered in our study. To analyze 
the effects of disorder and donor acceptor interfaces in 
more detail, Monte-Carlo simulations are used which are 
not the scope of this study [H, [27|, H, [H, [H SH [H 13 . 
Albeit all the simplifications, the Onsager-Braun theory 
is a widely used model to describe the dissociation of 
CT-excitons in organic semiconductors. 



2. Drift- diffusion modeling 

In the semiconductor, the following drift-diffusion 
equations apply: 



dE{x) 
dx 

J e (x) 

dn(x) 
dt 



e r eo 



\p(x) - n(x)} 



e/j e n(x)E(x) + D{^jl,T) 



dn(x) 
dx 



(7) 
(8) 



1 dJ e {x) 
e dx 



r e ffr(x)p(x)n(x) +kdS(x).(9) 



The electrical field E(x) in equation ([7]) is calculated 
dependent on the position x. Where e stands for the ele- 
mentary charge, n(x), p(x) for the electron/hole density, 
£o for the permittivity of free space and e r for the relative 
permittivity 

Equation ([8]) is the current equation for electrons. The 
mobility fj, e is shown to be independent of the electrical 
field, temperature or density, which is a simplified as- 
sumption used in this paper. The diffusion coefficient is 
given by the Einstein relation: D = iLksTjq. 

Equation © is the continuity equation for electrons. 
This equation takes the creation, migration and recombi- 
nation of the charge carriers into account. r(x) stands for 
the Langevin recombination r(x) = (/i e + I j >p)j^~- The 
requirement for this model to be applicable is that the 
mean free path A < r c , where r c stands for the coulom- 
bic capture radius r c = 47re e £ofer which for an average 
organic solid (eo ~ 3 - 4) gives a coulombic radius of 
r c rs 14 — 19nm. This requirement is satisfied since the 
mean free path in this low-mobility disordered organic 
semiconductor is smaller than the coulombic capture ra- 
dius [34[ and thus this model is applicable. The Langevin 
recombination efficiency prefactor, r e ff, is a factor be- 
tween and 1 which has been introduced by Juska et 
al. [35l ] to take the reduced Langevin recombination into 
account. It will be discussed later in this paper. 
For equations and © exist analogous equations for 
holes instead of electrons. The above set of equations is 
spatially discretised in SETFOS 19] using the Scharfet- 
ter Gummel discretisation [361 ] and solved iteratively in 
time using the Gummel scheme [37j . 



3. Built-in voltage 

There is a debate on the nature of the open-circuit 
voltage V oc in OPVs and currently there are several al- 
ternative explanations. There is evidence, that the open- 
circuit voltage depends on the highest occupied molecu- 
lar orbit (HOMO) and lowest unoccupied molecular orbit 
(LUMO) difference of the donor and acceptor molecules 
at the interface [38[ (3{|. Recently, this mechanism has 
been looked at in more detail and a linear dependence of 
the V oc on the energy of the charge-transfer absorption 
has been found [40(. The open-circuit voltage also de- 
pends on the work-function of the electrodes [4l[ . This 
influence seems to be more pronounced in bulk hetero- 
junction than in planar heterojunction device setups [42| . 
Furthermore the open-circuit voltage was observed to be 
intensity and temperature dependent 43] . This complex 
behavior is challenging for the simulation of these devices, 
because there is not a single theory which is capable of 
describing all of those dependencies. The approach used 
here assumes that every device has a built-in voltage 
which is derived from the difference in the work functions 
of the two electrodes. This is then used to calculate the 
effectively applied voltage V e / / . The electric field distri- 
bution inside the device is obtained through integration 
of the Poisson equation (cf. Eq. (O) using the effective 
bias V e // as a constraint to determine the integration 
constant, 



Veff = V app i - Vu = I Edx, 
Jo 



(10) 



where Y app i is the experimentally applied voltage. 
Therefore, the electric field due to space charges is su- 
perimposed on the applied field. 



4- Charge extraction 

When charges reach the organic-metallic interface they 
are extracted from the device. The current at such an in- 
terface has been described by Scott et al. [44[ as a balance 
of injection and surface recombination. This model con- 
siders the barrier reduction at an organic-metal interface 
due to the electric field and the image charge potential 
and calculates the net injection current. 



5. Validation of the simulator 

The model presented has been validated by varying the 
input parameters and analyzing their influence on the 
output. We have observed, that the open-circuit voltage 
increases slightly with illumination intensity. It also de- 
pends on the mobility of the charge carriers. Furthermore 
it is not equal to the built-in voltage which is assumed 
to be the difference between the work-functions of the 



Description 



Parameter Value Unit 



Workfunction PEDOT 
Workfunction Al 
HOMO P3HT [46] 
LUMO PCBM [46] 
Electron mobility [45] 
Hole mobility [45] 



5 eV 
4.31 eV 
5.3 eV 
3.7 eV 



5xl0~ 7 
lxlO~ 8 



2 /Vs 
2 /Vs 



TABLE I: Input parameters for the thickness dependent sim- 
ulation of figure [4] 



electrodes in this model. The fill factor is influenced by 
recombination losses and thus the mobility of the charge 
carriers and the Langevin recombination efficiency r e //. 
The short-circuit current depends linearly on the illumi- 
nation intensity until recombination losses take over and 
the increase becomes sub-linear. The results of this vali- 
dation correspond to experimental observations. 



III. ESTIMATION OF THE DISSOCIATION 
RATE 
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FIG. 4: Thickness dependence of the short circuit current 
Jsc with varying Langevin recombination efficiency r e ff. A 
purely optical simulation is also shown. The simulations are 
compared with measurements (■) by Gilot et al. [20( 



After describing the coupled opto-electronic nu- 
merical device model in the last section, it will be 
applied in this section to simulate an OPV device 
and extract unknown CT-exciton dissociation input 
parameters. The layer structure is shown in figure 
[TJ The indium tin oxide (1TO) layer has a thick- 
ness of 130 nm, it is followed by a 50 nm thick 
Poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) 
(PEDOT:PSS) layer. The active layer is modeled as a 
1:1 weight ratio P3HT:PCBM blend which is covered by 
100 nm Aluminum. The parameters for the electrical 
simulations which are taken from literature are shown in 
table |T1 The two mobilities are taken from Mihailetchi 
et al. [45| , who have measured the constant mobilities of 
electrons and holes in a P3HT:PCBM bulk heterojunc- 
tion solar cell depending on the annealing temperature. 
The measurements used here are for samples which have 
been annealed at 140° C for 5 minutes and the mobilities 
have been chosen accordingly. 

The unknown parameters are the decay rate kf, the 
pair separation distance a and the photon-to-CT-exciton 
conversion efficiency g e f / which have been introduced in 
equation ([5]) and ([6]) respectively. These parameters can 
be determined by applying the device model to thick- 
ness dependent measurements of the short-circuit cur- 
rent. Some reports in literature have used a simplified 
OPV model without CT-exciton dynamics (cf. Juska et 
al. [IBj]) and have found reduced recombination rates. 
Similarly our first analysis step assumes that absorbed 
photons directly generate free electron hole pairs and that 
recombining charges are lost and not fed into the conti- 
nuity equation for CT-excitons ©. This reduced model 
thus consists of equations ([7]), ([5]) and the modified ver- 
sion of equation © which reads as follows 



~df = e~ ±c r eff r(x)p(x)n(x) + g eff G opt (($!,) 

where g e ff stands for the photon-to-electron conver- 
sion efficiency. 

Figure [4] shows a purely optical simulation (top curve) 
to calculate the maximum achievable photocurrent 
depending on the device thickness. The other six curves 
are for coupled opto-electronic simulations where the 
Langevin recombination efficiency r e ff has been varied 
between 1 and 0.01. These simulation are compared 
with measurements by Gilot et al. [20j ■ For devices with 
an active layer thickness below 50 nm electrical losses 
can be neglected. Therefore the photon-to-CT-exciton 
conversion efficiency g e // is fitted to this part of the 
curve. A value of g e ff— 0.66 has been determined which 
is consistent with the analysis by Gilot et al. on the 
same thickness dependent data. For thicker devices elec- 
trical losses play an important role. Due to measurement 
uncertainties it is hard to give a final value, but figure 0] 
suggests that the Langevin recombination efficiency r e t / 
in the simplified model is 10% or lower. In the full model 
there is an equilibrium between CT-exciton formation 
and dissociation. We translate the observation of the 
reduced r e ff in the simple model to the dissociation 
probability in the full model, since the rate r x r e ff 
(simple model) corresponds to r x (1 — P) (full model) 
(cf. with equations © © )• Thus P must be 90% or 
higher. 

With this information it is possible to determine 
the unknown CT-exciton parameters by analyzing the 
Onsager-Braun model of CT-exciton dissociation. The 
initial pair separation distance a can be determined un- 
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FIG. 5: Dissociation probability according to the Onsager- 
Braun theory depending on the electrical field for several ini- 
tial pair separation distances a. The decay rate has been kept 
fixed at kf = 1 X 10 5 s" 1 . 



der the assumption that kf is set to 1 x 10 5 s _1 . The 
current density in figure |4] has been measured and cal- 
culated under short-circuit current conditions. From the 
open-circuit voltage which is about 0.6 V and the device 
thickness the internal electric field is estimated to be in 
the range of lxlO 6 to 2xl0 7 V/m for thick and thin 
devices respectively. Figure [5] shows the desired dissocia- 
tion efficiency (over 90%) and the internal electrical field 
range together with different initial pair separation dis- 
tances which have been varied between 1.5 nm and 1 nm. 
From this plot alone the only condition for a is that it has 
to be larger than 1.2 nm. The best fit value for the pair 
separation distance a has been chosen to be 1.285 nm, 
by comparing experimental current-voltage curves with 
simulated curves for an active layer thickness of 70 nm 
(cf. Fig. [6]). The thickness of 70 nm has been chosen 
because it is close to the first peak of the short-circuit 
current. Now, all the parameters for the full model are 
determined and summarized in table U and [TTJ They are 
used in the following simulations where the full model 
including CT-exciton dissociation is considered. 

Note that in the full model there is no need to reduce 
the recombination efficiency. Even though the original 
recombination rate is used (i.e. r e // = 1), there are small 
recombination losses due to the competing dissociation. 
Also note that the value of g e f / which was determined 
in the simple model is also applicable in the full model 
since with either model the parameter extraction is done 
by dividing the calculated theoretical maximum of the 
photocurrent by the measured one, preferably for thin 
devices were transport losses can be neglected. 



FIG. 6: Comparison of a simulated and experimentally mea- 
sured current- voltage curve for an active layer thickness of 70 
nm. The parameters used for this simulation are summed up 
in table 1 and 2. 



Parameter 



Symbol Value Unit 



Recombination eff. r eff 

Decay rate kf 
Pair separation a 

Optical generation eff. g e f / 



1 

1x10 s s _1 
1.285 nm 
0.66 



TABLE II: Derived simulation parameters from figures |4j [5] 
and [6] 



IV. SENSITIVITY ANALYSIS 

In the previous section a set of parameters has been 
found (cf. table Q] and |TT] ) such that the simulation fits 
the measurements. In this section the sensitivity of the 
model on the input parameters is investigated. This sen- 
sitivity analysis gives insight into simulation as well as 
experiment: First we analyze and study which param- 
eters influence the outcome of the simulation, secondly 
the sensitivity analysis also helps to identify experimen- 
tal setups where the influence of different parameters can 
be separated. 

The sensitivity analysis has been conducted as follows: 
Firstly, a reference has been simulated with the parame- 
ters indicated in table HI and ITT1 Secondly, five input pa- 
rameters have been chosen, Langevin recombination ef- 
ficiency r e ff, the electron/hole mobility /x„//x p , the pair 
separation distance a and the decay rate kf. The first 
three parameters (r e //, fJ-n, Mp) have been varied ± 10% 
and because the influence of the CT-exciton parameters 
(a, kf) is larger, these two parameters have been varied 
by 3% only. Therefore in this step ten simulations have 
been done (five parameters with variation in two direc- 
tions each). As a last step, the reference curve has then 
been subtracted from the variated curves. This difference 
between the reference curve and the varied one has then 
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FIG. 7: Sensitivity of the current density for a thickness 
dependent simulation. 



FIG. 8: Sensitivity of the current density in for a current- 
voltage curve. 



been plotted for each parameter. This procedure reveals 
possible orthogonalities among the parameters. 



A. Thickness dependent sensitivity 

First of all the thickness dependence of the short- 
circuit current has been analyzed, the results are shown 
in figure [7] The measurements are also shown, and agree 
the simulation for active layer thicknesses up to 150 nm. 
The two measurements around 230 nm which are shown 
in figure [5] are lying outside the axis scale in figure [7] and 
can not be modeled with the parameter set used in this 
study. This suggests that the different spin-speed and 
solvent concentrations that have been used to achieve 
the higher active layer thicknesses 20] also changed 
electrical properties of the active layer and therefore can 
not be modeled with a single set of parameters. 

As is shown in figure [71 the influence of the electric 
input parameters are only relevant for an active layer 
thickness larger than about 50 nm, which has already 
been seen for the Langevin recombination efficiency in 
figure[H The pair separation a has the highest influence 
on the short-circuit current although it has been varied 
by only 3%. The influence of all the parameters grows 
with increasing active layer thickness. 
Further analysis, not shown here, reveals a strong linear 
dependence between the three parameters pair separa- 
tion distance a, decay rate kf and electron mobility jj, n . 
These three parameters are all involved in the dissocia- 
tion of CT-excitons (cf. Eq. ^j). The thickness variation 
influences the electric field inside the device. Therefore 
the electric field dependence of the CT-exciton dissoci- 
ation is probed. The three parameters which show the 
same thickness dependence are all pre-factors for the field 
dependence as shown in Eq. In these equations the 
hole mobility is added to the electron mobility to deter- 



mine the dissociation rate kj, but in the device investi- 
gated here the electron mobility is higher by a factor of 
50 than the hole mobility and therefore dominates the 
dissociation. The hole mobility fi p and the Langevin re- 
combination efficiency r e / / both have a unique thickness 
dependence, which is slightly different from each other. 



B. Current-voltage curve sensitivity 

Figure [8] shows the sensitivity analysis for a current- 
voltage curve with an active layer thickness of 70 nm. 
The current-voltage curve is the same one as in figure 
[5] The sensitivity analysis clearly demonstrates that all 
parameters have the highest influence on the current- 
voltage curve in the fourth quadrant and close to the 
open-circuit voltage V oc . Again the pair separation dis- 
tance a has the strongest influence on the output. The 
hole mobility clearly has a different influence on the out- 
put than all the other parameters. Especially for applied 
voltages which are higher than V oc the hole mobility has 
the highest influence on the current. This suggests that 
the hole mobility is the limiting factor for the current 
(photocurrent plus injected current) drawn from the de- 
vice at higher voltages. 

Further analysis (not shown) reveals that for the 
current-voltage curve under consideration there are only 
two clearly linearly independent sets of parameters. The 
first set comprises the pair separation a, the decay rate 
kf, the electron mobility \i n and the Langevin recombi- 
nation efficiency r e f / , the second one is the hole mobility 
fi p . Although the Langevin recombination has a slightly 
different voltage behavior than the other three parame- 
ters in this group it is very similar. This is interesting 
because in the sensitivity analysis of the thickness depen- 
dent short-circuit current shown in figure [7] the influence 
of the Langevin recombination efficiency was more simi- 
lar to the one of the hole mobility. 
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FIG. 9: Normalized comparison between measured and sim- 
ulated J-V curves for different active layer thicknesses. 



FIG. 10: Sensitivity of the short-circuit current during turn- 
on. 



For all the measured short-circuit current points in fig- 
ure [4] a current- voltage curve was measured as well. An 
analysis of the influence of the device thickness on the 
current-voltage curve is performed for three representa- 
tive thicknesses which is shown in figure [5J The short- 
circuit current as well as the open-circuit voltage depend 
on the amount of the absorbed photons. Therefore to 
compare the current-voltage curves for different thick- 
nesses the short-circuit current as well as the open-circuit 
voltage need to be normalized to exclude optical effects. 
The comparison in figure [5] shows, that the shape of 
the experimental current-voltage curve does not change 
for different thicknesses of the active layers. In contrast 
the simulated current-voltage curves show a clear depen- 
dence of the shape on the device thickness. This indi- 
cates that in the fourth quadrant the influence of the 
electric field on dissociation efficiency and thus recombi- 
nation losses is overestimated in the device model pre- 
sented here. A possible explanation for this discrepancy 
might be the nature of the Onsager -Braun dissociation 
mechanism, which does not consider energetic disorder. 
Though much used, the Onsager-Braun model is not nec- 
essarily the best physical description of charge separation 
and more refined models have been developed. For in- 
stance Wojcik et al. [47| have derived an extended model 
based on the Onsager theory which leads to a weaker field 
dependence. 



C. Transient current sensitivity 

The sensitivity analysis in thickness dependence and 
current-voltage curves have shown, that several input 
parameters have a similar influence on the output and 
thus those parameters can not be determined from these 
measurements. In this part transient simulations are dis- 
cussed and their sensitivity on parameters is analyzed. 



Figure [10] shows the simulated transient of the short- 
circuit current when the light is turned on at t=0, for 
a thickness of the active layer of 70 nm. Again the pa- 
rameter set shown in table |T] and |TT] has been used. The 
sensitivity clearly shows that in the turn-on dynamics 
the Langevin recombination efficiency r e // and the de- 
cay rate kt do not play an important role, therefore the 
turn-on dynamics is solely defined by the pair separation 
distance a and the two mobilities. The CT-exciton decay 
rate and thus the recombination of free charge carriers 
play no role. The electron mobility determines the first 
part of the turn-on and the sensitivity reaches its maxi- 
mum at 0.15 /xs, whereas the hole mobility has its max- 
imum at 0.35 /is. Further analysis (not shown) reveals 
that the electron mobility /i n has exactly the same influ- 
ence on the turn-on behavior as does the pair separation 
a. This indicates that the turn-on behavior is limited 
by the time it takes to dissociate the CT-exciton pair. 
When the CT-exciton parameters are chosen such that 
the average lifetime of the CT-exciton is much smaller 
than the time it takes to extract the charge carriers, the 
turn-on behavior is solely defined by the two mobilities. 

Figure [IT] shows the transient of the short-circuit cur- 
rent when the light is turned off at t=0. The shape of 
the transient turn-off curve shown in figure [TT] is clearly 
different from the turn-on curve shown in figure 1101 The 
first one has a clearly visible kink at t = 0.015 /is whereas 
the latter curve is more smooth without any visible kink. 
The sensitivity analysis of the turn-off behavior shows 
that the curve is solely dependent on the two mobilities. 
This means that CT-exciton and recombination param- 
eters play no role and thus the recombination and sub- 
sequent CT-exciton dissociation play no role. The initial 
drop is given by the mobility of the electrons pb n whose 
influence peaks at t=0.01 us whereas the second, much 
slower drop is given by the hole mobility /z p whose influ- 
ence peaks at t=0.45 fis. This difference is in the same 
range as the difference in the mobilities of the two charge 
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FIG. 11: Sensitivity of the short-circuit current during turn- 
off. 

carriers which is — = 50. 

This last sensitivity analysis shows, that by measuring 
the short-circuit current during turn-off, the two mobility 
parameters can be estimated. This observation is similar 
to the findings by Pinner et al. 48] for electrolumines- 
cence transients in organic light-emitting diodes where 
the turn-off regime is governed by the built-in field as 
well. The time it takes to extract all charge carriers and 
thus the mobilities can be calculated under the assump- 
tion that the active layer is homogeneously filled with 
charge carriers and that the internal electric field is given 
by the difference on the applied and the built-in voltage 
divided by the active layer thickness E = ^ v ~Y bi ' , We 
thus obtain 



M = TTT ttt-- W 

(Vappl - VbijT 

The parameter t is the time it takes for the current to 
drop down to 5% of the initial value. With this approach 
the electron as well as the hole mobility can be extracted 
separately. 

Figure [12] shows the turn-off behavior as a semi- 
logarithmic plot. The first steep drop of the current is 
associated with the electrons and thus a line is fitted 
to it. Point A' is defined as the point along the fitted 
curve were the initial current at point A is reduced to 
5%. The same is done for the second drop which is as- 
sociated with the holes. Again B' is defined as the point 
where the current dropped down to 5% with B as the ref- 
erence point. This gives for t„=0.037 [is and r p =0.56 /is. 
With V app ;-Vki=0.69 V (cf. table |J) and an active layer 
thickness 70 nm the following mobilities are extracted: 
fi n = 1.9xl0~ 7 m 2 /Vs and fi p = 1.3xl0~ 8 m 2 /Vs. The 
extracted values are in the same order of magnitude as 
the input values given in table [I] The mobility of the fast 
charge carriers is underestimated and the mobility of the 



FIG. 12: Semi-logarithmic plot of the turn-off behavior shown 
in figure ITTI This plot can be used to extract the mobilities 
of electrons and holes. 



slow charge carriers is slightly overestimated. This dis- 
crepancies stem from the assumption of a homogeneous 
electric field inside the device and thus the independent 
treatment of electrons and holes which is a simplification. 



V. CONCLUSION 

From the measured thickness-dependent short-circuit 
current for thin active layers and the optical model we 
estimated a photon to CT-exciton conversion efficiency 
g e ff of 66%. Using this dataset the coupled opto- 
electronic simulation suggests a lower limit for the CT- 
exciton dissociation efficiency of 90%. Adding the mea- 
sured current-voltage curve to the numerical analysis and 
assuming that kf is equal to 1x10 s s _1 we find dissoci- 
ation parameters that are consistent with the estimated 
90% dissociation efficiency. From a comparison of several 
current- voltage curves at different thicknesses there is ev- 
idence that the influence of the electric field on the CT- 
exciton dissociation process is overestimated using the 
Onsagcr-Braun model. Further investigations have to be 
carried out to clarify this issue. The sensitivity analysis 
conducted demonstrates that the influence of the two ex- 
citon parameters and the electron mobility are linearly 
dependent in the current- voltage curve and photocurrent 
thickness scaling. A sensitivity analysis of transient pro- 
cesses indicates that the input parameters can be sep- 
arated. During turn-on, the CT-exciton decay rate kf 
and thus the charge carrier recombination play no role. 
Even more insight can be gained by looking at the turn- 
off behavior. In this regime the transient current density 
is solely defined by the mobilities of the charge carriers. 
This suggests that the charge carrier mobilities can be 
extracted. The influence of charge traps was not consid- 
ered in this study but will be discussed elsewhere. 
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